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I .   INTRODUCTION 

Today,  naval  architects  are  greatly  interested  in  the  design  of 
super  tankers .   However,  many  difficulties  arise  in  design  when  the 
strength  of  the  super  tanker  is  considered  from  an  economic  point  of 
view.   The  transverse  strength  problem  is  apparently  one  of  these  diffi- 
culties . 

There  are  several  ways  in  analysing  the  transverse  strength  prob- 
lem.  In  each  method,  the  two-dimensional  approach  is  simpler,  but  the 
three-dimensional  approach  is  so  complicated  that  hardly  no  one  wishes 
to  use  it.   Therefore,  some  approximate  is  devised  and  sufficient  results 
are  achieved.   But  still  some  architects  wish  to  know  the  three- 
dimensional  effect  on  transverse  strength  analysis. 

Some  prominent  naval  architects  suggested  solving  the  three- 
dimensional  strength  calculation  but  most  of  them  did  not  apply  numerical 
computations  to  existing  complicated  ships,  for  even  though  the  computor 
is  available,  it  is  still  no  easy  task.   So  many  variables  exist  it  is 
very  tedious  work  to  put  them  into  computors . 

In  this  paper  we  are  going  to  consider  these  tedious  problems  by 
using  FORTRAN  7090. 

As  previously  mentioned,  there  are  many  ways  to  solve  the  three- 
dimensional  transverse  strength  problem  and  each  one  ha.s  its  own  merits, 
although  they  will  give  almost  the  same  level  of  accuracy. 

Since  some  approaches  are  so  complicated  and  give  the  same  level 
of  accuracy  as  more  simple  approaches,  they  become  useless.   Therefore, 
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we  are  going  to  deal  with  a  well-known  and  more  simple  approach--the 
slope-deflection  method--to  attack  these  three-dimensional  calculations 
and  examine  its  effect. 

In  large  tankers  especially,  each  span  has  such  deep  depth- 
breadth  ratios  that  we  suppose  shear  force  contributes  more  than  usual, 
therefore,  we  take  into  account  shear  force  effect  in  slope  deflection 
equation. 

Eventually,  we  will  check  our  computed  results  by  means  of 
slope-deflection  method  with  one  of  the  fixed  point  methods  which  was 
given  in  our  class  by  Professor  Schade.   This  fixed  point  method  differs 
in  assumption  from  the  slope-deflection  method,  but  we  can  make  this 
fixed  point  assumption  comparable  with  our  slope-deflection  assumption 
by  restricting  our  loading  conditions.   Although  we  are  dealing  with  the 
fixed  point  method,  it  will  not  give  us  shear  force  effect  and  three- 
dimensional  effect.  We  can  see,  therefore,  only  the  differences  between 
those  two:   the  slope-deflection  method  with  shear  force  effect  and  the 
fixed  method  without.   They  will  lead  us  to  some  interesting  conclusions, 


II.   BASIC  SLOPE -DEFLECT ION  EQUATION  CONSIDERED  WITH  SHEAR  EFFECT 
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(2)   Simple  beam  with  bracket  at  both  ends 
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Using  all  of  these  results  and  the  result  of  (l),  we  can  easily  get 
the  following  basic  formula  using  the  following  notation. 
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and  the   constant  terms      C^    ,    Cnffl   ,   D^    ,    Dnm  are: 
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For  the  purpose  of  our  calculations,  we  only  need  two  kinds  of 
loading  spans.   They  are  shown  on  Table  I  and  Table  II. 
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III.   SETTING  UP  SIMULTANEOUS  EQUATION 

A.  Assumption 

Once  the  designer  has  made  the  decision  to  consider  the  trans- 
verse frame  of  a  ship,  he  is  then  faced  with  the  task  of  deciding  on 
simplifying  assumptions.   The  main  objective  of  these  assumptions  is  to 
create  a  fictitious  structure  and  condition  of  loading,  which  will  re- 
duce the  amount  of  labor  in  performing  the  structural  analysis  without 
compromising  too  severely  the  validity  of  the  numerical  answers. 

(1)  One  of  the  factors  a  designer  must  consider  is  the  degree  of 
dimensional  similitude  he  can  retain  in  the  structure  he  selects  for  his 
computational  work.   For  example,  he  must  decide  how  much  of  the  shell 
plating  should  be  included  in  the  analysis  of  a  transverse  frame  as  ef- 
fective material.  We  are  going  to  use  Effective  Breadth  Conception  in 
Ship  Structure  Design,  published  in  195^-  by  Professor  Schade. 

There  are  several  other  assumptions  to  be  applied  in  this  paper. 

(2)  The  relative  angle  initially  formed  by  the  tangents  to  the 
elastic  curve  at  a  point  where  two  or  more  members  are  joined  are 
assumed  to  remain  unchanged  after  loading;  that  is,  if  two  members  are 
originally  joined  at  right  angles,  the  angle  formed  by  the  tangents  to 
the  elastic  curve  after  the  loads  are  applied  will  remain  90°- 

(3)  Joints  of  1-1  (or  9-2)   and  k-2    (or  13-2)  are  fixed  in  space 
but  free  to  rotate  in  three-dimensional  analysis  and  additionally  9-l 
(or  l^--2),  15-2  (or  13-1),  14-1  and  15 -1  are  fixed  in  space  in  two- 
dimensional  analysis  (see  Fig.  2-2). 

B.  Simultaneous  Equation 

Before  setting  up  simultaneous  equations  we  must  investigate  the 
number  of  unknown  variables.  From  the  basic  slope  equation  we  can  take 


up  the  following  variables  as  unknown  variables . 

1.  91(1)  =  92(9)  =  X(l) 

2.  91(9)  =  ©1(5)  =  92(14)  =  X(2) 

3.  92(1)  =  92(10)  =  91(2)  =  X(3) 

4.  92(5)  -  91(10)  =  91(6)  =  X(4) 

5.  92(2)  =  92(11)  =  91(3)  =  X(5) 

6.  92(6)  =  91(11)  =  91(7)  =  X(6) 

7.  92(3)  =  92(12)  =  91(4)  =  X(7) 

8.  92(7)  =  91(12)  =  91(8)  =  X(8) 

9.  92(4)  =  92(13)  =  X(9) 

10.  92(8)  =  91(13)  =  92(15)  -  X(10) 

11.  jrf(l)  =  i^(5)  =  X(ll) 

12.  t(2)   =^(6)  =  X(12) 

13-  jrf(3)  =  Ml)   =  x(!3) 

14.  ^(4)  =  jrf(8)  =  X(l4) 

15.  H9)   =  ^(10)  =  ^(11)  =  ^(12)  =  ^(13)  =  X(15) 

16.  <j>(±k)   =  X(l6) 

17.  ^(15)  =  X(17) 


Now  we  have  17  unknown  variables  in  one  transverse  frame. 
Therefore  we  have  to  have  at  least  17  simultaneous  equations  for  each 
frame . 

1.  Ml(l)  +  M2(9)  =  0 

2.  Ml(9)  +  Ml(5)  +  M2(l4)  =  0 
3-   M2(l)  +  M2(10)  +  Ml(2)  =  0 

4.  M2(5)  +  M1(10)  +  ML(6)  =  0 

5.  M2(2)  +  M2(ll)  +  Ml(3)  =  0 

6.  M2(6)  +  ML(ll)  +  Ml(7)  =  0 
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7-  M2(3)  +  M2(12)  +  ML(k)  =  0 
8.  M2(7)  +  Ml(l2)  +  ML(8)  =  0 
9-   M2(4)  +  M2(13)  =  0 

10.  M2(8)  +  Ml(l3)  +  M2(15)  =  0 

11.  01(6)  +  Ql(2)  -  Q2(l)  -  Q2(5)  =  0 

12.  Ql(7)  +  01(3)  -  02(2)  -  Q2(6)  =  0 
13-  01(8)  +  Ql(10  -  Q2(7)  -  02(3)  =  0 

Ik.      X(ll)  x  SPAN(l)  +  X(12)  x  SPM(2)  +  X(l3)  x  SPAN(3)  +  X(l^) 

x  SPAN(1+)  =  0 
For  the  three  dimensional  analysis  we  can  set  up  the  following 
three  more  equations . 

15.  YB  =  X(15)  x  SPAN(13)  +  X(17)  x  SPAN(l5) 

16.  YD  =  X(15)  x  SPAN(9)  +  X(l6)  x  SPAN(l4) 

17.  YL  =  X(15)  x  SPAN(13) 

Here  Y-n  ,  Y-q  and  Yx  are  deflection  of  bottom  and  deck  longitudi- 
nals and  longitudinal  BHD.  We  suppose  that  longitudinals  are  only  acted 
by  the  concentrated  load  at  joints  with  transverse  frames.   Then  acting 
load  P-n  ;  Prv  and  Pj  are  expressed  as  follows: 

PB  =  201(15) 

PD  =  2Q1(U) 

PL  =  01(9)  +  01(10)  +  01(11)  +  01(12)  +  01(13)  -  02(15)  -  02(1*0 
For  the  two  dimensional  analysis  we  need  only  to  solve  Ik   simultaneous 
equations  with  unknown  variables  Ik   by  setting  simply  X(l5)>  X(l6)  and 
X(l7)  to  zero  in  the  above  equations. 

In  this  paper,  for  the  practical  numerical  computation,  we  will 
take  an  existing  Japanese  built  Tanker  as  an  example.   She  has  ^7^000 
tons  Dead  Weight  and  one  of  the  tank's  length  is  503- 936" ]    and  height 
60l+.329"i  and  its  breadth  is  II88.97V.   One  of  the  web  frame** 
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scantling  is  shown  in  Fig.  1  and  three  web  frames  placed  in  a  length  of 
tank. 
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IV.   PRELIMINARY  CALCULATION 

A.   Section  modulus  calculation 

From  the  effective  "breadth  conception  curve  CASE  (III),  we  can 
easily  find  out  P   :  Effectiveness  ratio.   From  these  data  we  can  calcu- 
late moment  of  inertia  and  neutral  axis. 

The  results  are  as  follows: 


SPAN  No. 

P 

MOMENT   OF   INERTIA 
OF  THE  SECTION 

RATIO  OF 
INERTIA 

DISTANCE  FROM  N.A. 
TO  FREE  FLANGE 

1 

0.570 

30828.851 

0.913 

47.21+0 

2 

0.350 

25839.827 

O.79I+ 

42.339 

3 

0.350 

25839.827 

O.79I+ 

42.339 

4 

0.44-0 

1+5197.227 

0.822 

^1. 335 

5 

0.635 

23601+.205 

O.89I 

42.189 

6 

0.410 

20716.733 

O.78O 

38.867 

7 

0.1+10 

21011+.376 

O.78I+ 

39-184 

8 

0.485 

6081+7.290 

0.814 

33-413 

9 

O.870 

22181.710 

O.985 

44.287 

10 

1.000 

3083.156 

1.000 

16.259 

11 

1.000 

3656.705 

1.000 

16.377 

12 

1.000 

5282.3^7 

1.000 

16.748 

13 

0.870 

8I+7I+1 .  3I+6 

0.971+ 

55.841 

Ik 

0.735 

21+303.928 

0.963 

43.268 

15 

0.735 

60303.312 

0.91+5 

57-421 

Table  3 
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In  inches 


SPAN 
No. 

SPAN  LENGTH 
L 

t 

h 

l2 

1 

198.454 

122 .  859 

58.367 

17.228 

2 

125.984 

9O.967 

17.372 

17.645 

3 

125.984 

89.928 

17.645 

18.411 

4 

137.987 

40.587 

17-843 

79.554 

5 

198.454 

120.428 

60.574 

17.452 

6 

125.984 

90.536 

17.596 

17 . 892 

7 

125.984 

89 .  456 

17.857 

18.671 

8 

137-984 

39.287 

17.955 

80.745 

9 

301 . 006 

176.205 

61.037 

63.764 

10 

301 . 006 

196.093 

52.051 

52 . 862 

ll 

301 . 006 

196.840 

51.670 

52 . 496 

12 

301.006 

198.488 

50.850 

51.668 

13 

301.006 

190.620 

54.579 

55.807 

14 

280.539 

233.377 

11 . 062 

36.IOO 

15 

280.539 

236.469 

12.451 

31.619 

Table  4 
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Fig.    2-1     Designation  of  span  number 
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Fig.  2-2  End  point  designation 
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V.   SOLVING  A  SIMULTANEOUS  EQUATION  USING  COMPUTOR  7O9O 

Before  solving  a  simultaneous  equation  we  decide  to  solve  two 
different  loading  conditions:   (l)  favorable  loading  condition  (see  Fig 
3),  (2)  most  critical  loading  condition  (see  Fig.  h) . 

For  the  girder  deflection  we  must  also  think  about  shear  defor- 
mation because  of  having  relatively  short  span  length  and  deep  depth. 
For  this  purpose  we  can  find  deflection  equations  of  simple  beam  at  the 
condition  of  both  ends  fixed  and  concentrated  at  some  span  point. 

One  can  use  the  next  simple  equation  for  the  concentrated  load 
at  the  center  of  the  beam. 


W^  =  -  -1-.  ^L?  (1  +  C^ 

kQ   E  ib 


P     192  By     P' 


C_  = 


b2  G  A> 


^b 

When  the  concentrated  load  is  acting  at  the  point  of  l/k   span 
length,  deflection  can  be  expressed  as  follows: 

W   =  -  J-  %   *!*  (1  +  2C  ) 
P     512  8  E  lh       Py 


Where  W  . . .Deflection  at  the  point  of  load  applied 
A^.  .  .Web  area  of  b  beam 
P  .  .  .Concentrated  load 
For  the  computed  deflection,  see  curves  ^— 1,  k-2   and  k-3> 
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One  of  the  computed  results  is  tabulated  below  as  Table  5 


^"""^^Angle 

Deflection  angle 

in  Radian 

X     ""\*^ 

x  10"^ 

X(l) 

0.475742 

X(2) 

-0.413502 

X(3) 

-0.226263 

x(4) 

0.128071 

X(5) 

0.171848 

x(6) 

-O.369586 

x(7) 

-O.501310 

x(8) 

O.550122 

x(9) 

0.647280 

x(io) 

-O.336503 

x(u) 

0.028208 

X(12) 

-0.080296 

X(13) 

-0.070333 

X(l4) 

0.096957 

Table  5 •   Acquired  unknown  variables  X  at  the 
forward  frame  (load  l)  in  two 
dimensional  analysis. 


From  these  values  we  can  get  moment  at  edge  1  and  2  (left  and 
right  side  or  upper  and  lower  side)  for  every  span.   You  can  see  one  of 
the  results  in  Table  6. 


2'3- 


SPAN  No. 

MOMENT  1  x  103 

MOMENT  2  x  1C>3 

1 

-o. 21+2163 

0.914909 

2 

-O.89I+309 

1.194515 

3 

-1.177446 

O.960957 

k 

-0.933445 

10.523114 

5 

0.388150 

-0.975332 

6 

0.936319 

-I.I5802O 

7 

1.201756 

-0.847029 

8 

0.799642 

-II.777653 

9 

-0.059339 

0.242128 

10 

O.039037 

-0.020693 

n 

-0.043852 

-0.170775 

12 

0.0^7344 

-0.027552 

13 

12.349345 

-10.523304 

14 

-0.169027 

-0.328803 

15 

-0.283449 

-0.571738 

Table  6.   Moment  of  the  forward  frame 

From  these  end  moments  we  can  get  very  easily  a  middle  portion 
moment  by  simply  superimposing  loading  to  each  span.   This  work  can  be 
done  very  easily  if  we  use  the  Computor  for  every  span.   The  results 
are  plotted  in  figures . 


VI.   FIXED  POINT  METHOD 
A.   General 

For  the  purpose  of  checking  above  results  and  investigating 
shear  effect,  we  will  apply  here  numerical  computation  by  the  method  of 
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fixed  point.  The  fixed  point  method  is  originated  from  the  three  moment 
equation.   In  this  approach,  we  have  the  following  assumptions: 

(1)  All  nodal  points  are  fixed  in  space. 

(2)  No  shear  force  is  considered. 

We  will  not  discuss  here  in  detail  about  this  derivation  (see  Appendix), 
but  will  apply  this  method  in  practical  calculation  using  the  same  load- 
ing conditions  to  the  same  scantling  frame  ship. 
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Then  the  moment  at  the  ends  of  the  span  can  be  expressed  by  the  follow- 
ing formulas . 
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£'  :  is  the  span  length 
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And  Re  and  Le  can  be  found  out  in  some  handbooks 


That  is 


L  e  =  V 

Re=  Ao(1"^) 

A  is  the  area  below  the  simple  moment  curve. 

and  x  is  the  distance  from  the  left  side  to  the  center 


of  the  area  AQ  . 
For  our  computation  purpose  we  can  find  out  the  following  data  from  the 
handbook. 
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And  Re  and  L£  can  be  found  out  in  some  handbooks 
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A  is  the  area  below  the  simple  moment  curve. 

and  x  is  the  distance  from  the  left  side  to  the  center 

of  the  area  AQ  . 

For  our  computation  purpose  we  can  find  out  the  following  data  from  the 
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If  once  we  get  moment  for  each  span,  then  we  can  compute  each  joint 
moment  using  stiffness  and  fixed  points  by  one-by-one  adding  process  as 
follows : 


But  in  doing  so,  we  must  be  very  careful  in  the  sign  for  some  opening 
span.   (Ref.  illustrative  calculation  of  moment  distribution  Table  12.) 
Now  here  is  the  practically  computed  results  listed. 


I.  k 

I. 

In 
^n 

■^n    ^15 

Unit  g       x  -=— ^~ 

c-n    £15 

1 

30828.851 

280.539 

2lk. 955 

2.1+80 

2 

8k7^1 .  710 

301.006 

281.528 

3-21+8 

3 

k5197-227 

137.987 

327.5^6 

3-779 

k 

608k7.290 

137.987 

kko.963 

5.O88 

5 

5282 . 3k7 

301 . 006 

17.5^8 

0.202 

6 

25839.827 

125.98k 

205.10k 

2 .  366 

7 

2101k. 376 

125.98k 

166.801 

1.92k 

8 

3656.705 

301.006 

I2.lk8 

o.iko 

9 

25839.827 

125.98k 

205.10k 

2 .  366 

10 

20716.733 

125.98k 

164. k39 

1.897 

n 

3083.156 

301 . 006 

10.2k2 

0.118 

12 

30828.851 

198. k5k 

l55-3k5 

1.792 

13 

2360k .205 

198. k5k 

118. 9k0 

1.372 

lk 

22181.710 

301.006 

73.691 

O.85O 

15 

2k303-928 

280.539 

86.632 

1.000 

Table  7 
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In  the  next  table  the  computed  values  of  A- n/Cn    ,    ^n/K  ,  ^/^ 
and  W  /k  are  shown . 


V^n 

6n/K 

S  It 

n/^n 

%/K 

1 

0.232 

1-055 

O.236 

1.022 

2 

O.219 

1-375 

0.224 

1.382 

3 

(0.250)0.248 

1.419 

9.256 

1.428 

4 

0.200 

2.159 

0.203 

2.164 

5 

0.240 

0.813 

0.245 

0.818 

6 

0.223 

1.107 

0.221 

1.104 

T 

(0.250)0.275 

1.154 

0.277 

1.189 

8 

0.240 

(3.816)3.021 

0.205 

2.920 

9 

0.327 

0 .  561 

0.262 

0.516 

10 

0.327 

0.077 

0.328 

0.077 

n 

O.327 

0.092 

0.328 

0.092 

12 

0.327 

0.135 

0.328 

0.133 

13 

0.260 

I.969 

0.222 

1.893 

14 

0.266 

0.610 

0.268 

0.612 

15 

0.284 

0.703 

0.285 

0.703 

Table  8 


On  the  next  page,   ///"   and  K/ A'   tables  are  also  prepared. 
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/>//>' 

Ay  A' 

1 

o .  308 

0 .  302 

2 

0,288 

0.208 

3 

0.3^ 

0.329 

4 

0.254 

0.250 

5 

0.324 

0.315 

6 

0.283 

0.287 

7 

0.383 

0 .  379 

8 

0.257 

0 .  315 

9 

0.355 

0.48.5 

10 

0.424 

O.485 

11 

0.424 

O.485 

12 

0.488 

0.485 

13 

O.285 

O.351 

14 

O.366 

0.362 

15 

O.398 

0,396 

Table  9-1 
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The  computed  end  moment  by  its  own  span  load  are  listed  in 
Table  10.   (We  restricted  our  computation  in  the  forward  frame  only.) 


SPAN  No. 
1 

LOAD 

I 

LOAD 

II 

MOMENT  1  x  103 

MOMENT  2  x  103 

MOMENT  1  x  103 

MOMENT  2  x  103 

0 . 42470 

0.52657 

0 . 42470 

0.52657 

2 

0.47568 

0.48144 

0.47568 

0 . 48144 

3 

0 . 81965 

O.79450 

O.81965 

O.79450 

4 

I.O3658 

1.03886 

I.O3658 

I.03886 

5 

-0.44029 

-O.54524 

0 

0 

6 

-0.46159 

-O.49807 

0 

0 

7 

-0.87103 

-O.89650 

0 

0 

8 

-O.95682 

-O.13608 

0 

0 

9 

0 

0 

0 

0 

10 

0 

0 

0 

0 

li 

0 

0 

0 

0 

12 

0 

0 

0 

0 

13 

-7.95302 

-5-59423 

-7.95302 

-5.59^23 

14 

0 

0 

0 

0 

15 

0 

0 

-6.82683 

-6.89544 

Table  10 


From  the  above  computed  end  moment  at  edge  1  and  2  we  can  dis- 
tribute moment  by  using  stiffness  and  fixed  points  f*    and  A^  •    (See 
Table  12.)   Then  we  can  get  the  following  final  end  moment. 
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Table  11.   Final  end  moment  in  tons-inch 


SPAN  No. 

LOAD 

1 

LOAD 

11 

MOMENT  1  x  103 

MOMENT  2  x  103 

MOMENT  1  x  103 

MOMENT  2  x  103 

1 

0.2968 

O.7627 

0.3396 

O.7616 

2 

0.7425 

1.1749 

0.7536 

1-1315 

3 

1.2091 

0.4132 

1.1643 

O.5683 

4 

0.3464 

6.2266 

0.5379 

5.6970 

5 

-0.3354 

-0.7174 

O.1312 

0.2012 

6 

-0.6946 

-1.3596 

0.2248 

-0.8502 

7 

-1.4136 

-0.1739 

-O.9191 

2.4365 

8 

-0.0980 

-7-3683 

2.6901 

-1.9521 

9 

0.1419 

O.2698 

-0.0935 

0.3396 

10 

-0.0228 

-0.0202 

-0.0236 

-0.0080 

11 

0.0535 

0 . 0342 

0 . 0689 

-0.0328 

12 

-0.0759 

-0.0668 

-0.2536 

-0.0304 

13 

-8.4853 

-6.2266 

-10.2938 

-5.6970 

14 

0 . 0662 

-0.1935 

0 . 0049 

0.0377 

15 

0.6203 

-1.1170 

-6.3110 

-8.3417 
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Table  12.   Illustrative  calculation  of  moment  distribution 
for  span  No .  1  load 
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-25.3 

63.6 

150.7 

-k2k.7 

-87.1 

27.1+ 

-12.3 

-5.1 

0.6 

22.3 
26.8 

-11.7 

26.5 

7.2 

3^-0 

3.8 

-7-7 

1.9 

-1-3 

-0.6 

-0.1 

3-1 

-9.0 

-7-1 

-7.1 
-1.9 

2.7 

-1.2 

-0.4 

-10.2 

2.3 

-1.2 

0.8 

0.2 

2.5 

-1.0 

2.5 

3-3 

2.5 

-0.8 

1.7 

-0.3 

1.0 

3-5 

-10.1 

0.3 

-0.5 

1.3 


3-0 


+ 
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(b) 

•+ 


\2k.7 


■526.5 

-500.0 


1^0  .0 

131.0 


-1+3-0 
40.5 


10.1 


(a) 


+ 
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Once  we  get  each  span's  end  moment  then  we  can  get  various 
point  moments  very  easily  by  superimposing  process.   And  these  numeri- 
cal computational  results  are  plotted  in  figures . 

Also  stress  curves  are  available  in  this  present  paper,  but  we 
did  not  compute  the  moments  and  stresses  for  all  the  frames  by  the 
fixed  point  method.  We  only  calculated  the  forward  frame  for  both 
loading  conditions . 

The  reason  is  this:   the  fixed  point  method  is,  of  course,  a 
two  dimensional  one  and  each  nodal  point  is  fixed  in  space;  therefore, 
a  little  difference  of  water  head  will  bring  us  a  very  insignificant 
computational  difference. 

VII.   DISCUSSION 

Most  of  the  naval  architects  are  anxious  to  learn  about  the 
three-dimensional  effects  on  the  transverse  frame  strength  of  the  ship. 
But  the  computation  is  rather  lengthy  and  tedious  to  work  out  so  the 
architects  estimate  the  effect  on  transverse  frame  strength  by  the 
simplified  method  which  satisfies  them.   In  this  paper  we  dealt  with 
this  lengthy  calculation  by  use  of  the  computor.   It  cannot  be  done 
without  this  machine  because  so  many  unknown  variables  are  interacted. 

Even  a  casual  glance  at  the  graph  will  show  that  there  are  not 
many  differences  between  the  two-dimensional  and  the  three-dimensional 
analysis . 

Especially  in  load  condition  |_,  the  results  are  nearly  the  same. 
Also  in  load  condition  2,  we  do  not  see  much  difference  between  these 
two . 

In  load  condition  1,  at  the  forward  frame  (Fig.  5 )  we  can  see 
also  almost  no  differences  between  the  fixed  point  method  and  our 
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present  method  in  the  upper  half  of  the  frame  that  is  span  No.  1,  2,  5, 
6,  9;  10;  11  and  Ik.      That  means  that  lightly  loaded  spans  do  not  have 
any  shear  effect  at  all,  but  loading  is  steadily  increased,  then  the 
discrepancy  is  noticeable. 

We  are  now  going  to  examine  these  graphs  individually  for  more 
detail. 

A.   Load  Condition  1 

(1)  This  loading  condition  is  ballasted  with  sea  water  in 
the  center  tank  and  the  side  shell  is  exhibited  to  the  hydrostatic  load. 
This  is  one  of  the  most  favorable  loading  conditions .   In  this  loading 
condition,  the  forward  frame  does  not  show  us  much  discrepancy  between 
two-dimensional  and  three-dimensional  analysis.   In  this  loading  condi- 
tion, junction  points  1-2  and  2-1;  2-2  and  3-1;  3-2  and  k-1;    also  5-2 
and  6-1;  6-2  and  7-1.5  and  7-2  and  8-1  remain  in  space  because  of  the 
symetric  loading  condition.   Therefore  the  fixed  point  method  and  our 
slope-deflection  method  can  be  comparable  in  quantities . 

Thus,  this  means  if  we  neglect  shear  effect,  we  will  lose  al- 
most U,000  Tons-inch  moment  at  the  span  edge  which  is  most  heavily 
loaded. 

(2)  Middle  frame:   We  can  see  a  little  difference  between 
these  two  methods,  but  not  too  much.   This  means  that  the  middle  frame 
is  the  most  effected  frame  among  the  three  frames  because  the  girders 
will  deflect  most  at  the  center. 

(3)  In  the  aft  frame,  we  can  see  almost  no  effect  of  the 
three-dimensional  effect  on  the  moment  curve.   Almost  the  same  as  the 
forward  frame . 
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B.   Load  Condition  2 

(1)  We  think  this  loading  condition  is  most  critical,  but 
still  two-dimensional  analysis  and  three-dimensional  analysis  have  the 
same  aspect.   They  don't  change  very  much,  but  here  are  several  inter- 
esting aspects  being  displayed.   That  is,  in  our  method  we  fixed  joint 
1-1  and  9_2  and  k-2   and  13-2  so  that  as  a  whole  side  frame  and  BHD  bent 
something  like  a  semi -fixed  beam  at  the  bottom  and  deck,  but  by  the 
fixed  point  method  nodal  points  are  fixed  in  space  so  that  they  have 
almost  the  same  aspect  as  in  load  condition  1.   But  from  this  curve  we 
can  see  also  a  big  difference  in  heavily  loaded  lower  spans. 

(2)  In  the  middle  frame  we  can  see  some  conceivable  dis- 
crepancies between  those  two  methods,  but  still  in  good  shape;  in 
three-dimensional  analysis,  however,  the  moment  at  point  15-1  and  9"2 
is  reversed  in  direction  because  of  the  restraints  by  girders  released 
at  15-1  and  13-1  points.   The  other  spans  have  almost  the  same  results, 
nevertheless  one  can  see  a  little  higher  moment  in  the  frames  9;  10;  Hj 
12  and  Ik. 

(3)  In  the  aft  frame,  the  aspect  is  almost  the  same  as  in 
the  forward  frame.  We  cannot  see  any  noticeable  discrepancies  between 
those  two  approaches. 

CONCLUSION 

Most  designers  of  super  tankers  are  worried  about  the  three- 
dimensional  effect  on  the  transverse  frame  strength.   But  from  our  com- 
putation we  cannot  find  large  discrepancies  between  two-dimensional  and 
three-dimensional  analysis.   However,  in  the  three-dimensional  analysis, 
calculation  is  more  detailed  and  tedious  than  in  the  two-dimensional. 
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For  our  present  discussion,  three-dimensional  computation  is 
more  academic  than  practical,  otherwise  we  would  not  use  any  simplified 
design  process  such  as  Professor  Schade's  Design  Curves  for  Cross 
Stiffened  Plating. 

Thus  our  conclusion  is:   forget  about  the  three-dimensional 
analysis,  unless  you  intend  to  use  any  simplified  design  process  pre- 
pared for  the  three-dimensional  analysis;  hut  take  into  account  the 
shear  effect  on  the  two-dimensional  analysis  on  super  tanker  designs. 


FIG.  5  FORWARD  FRAME  MOMENT  CURVE  (LOAD  1) 
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FIG.  6  MIDDLE  FRAME  MOMENT  CURVE  (LOAD  1) 
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FIG.  7  AFT  FRAME  MOMENT  CURVE  (LOAD  \) 
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FIG.  8  FORWARD  FRAME  MOMENT  CURVE  (LOAD  2) 
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FIG.  9  MIDDLE  FRAME  MOMENT  CURVE  (LOAD  2) 
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FIG. 10  AFT  FRAME  MOMENT  CURVE  (LOAD  2) 
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FIG.  11  FORWARD  FRAME  STRESS  CURVE  (LOAD  1) 
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FIG.  12  MIDDLE  FRAME  STRESS  CURVE  (LOAD  1) 
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FIG.  13  AFT  FRAME  STRESS  CURVE  (LOAD  1) 
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FIG*  14  FORWARD  FRAME  STRESS  CURVE  (LOAD  2) 
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FIG.  15  MIDDLE  FRAME  STRESS  CURVE  (LOAD  2) 
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FIG.  16  AFT  FRAME  STRESS  CURVE  (LOAD  2) 
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IX .      APPENDIX 

1 .      The  three  moment  equation 


i     pTTTTTTTTTT]  /<fTi  II  1  J,  111  IK      I    i~l  4  4  i  I I 

iff?  -,      t*V  "  TIT  _  7>V 


9~   = 


or 


7C 
tLC 


L 


L 


Then  by  using  the  second  moment  area  theorem  to  obtain 


tLC   and  ^C 


LC  "  Et 


Mxdx 


Therefore  simply 
1   1 


1       T        «  1 


- —  (Ax     +£lm-l    +±lm    -  l  ) 

L    EL        L  L        2      L  L  3      L        2      L  C    3      1/ 


LLEIL 


-  ^  4  (Vr  +  2  ^  3  **  +  2  ^  I  Lr) 


Simplifying  the  above  expression,  the  three  moment  equation  is: 


LR 


LR 
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2 .   Fixed  point  method 


By  three  moment  theory  (Ref .  Appendix  l)  first  three  support  (from  left) 
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v  xb    ^-c 
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cd 


=  0 


(where  AL  =  AR  =  0   ) 

Set    ty.  =  t-u   —   , 


Then 
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0  +  2Mbc  (£b  +  4)  +  Mcd  lc   =  0 
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By  geometry 
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Therefore  finally 
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Apply  three  moment  equation  to  span  C  and  D 
%ctc   +  2Mcd  (lc   +  ld)    +  Mde  ^d=  0 


using  relationship 
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In  general 

^n_   ^ 

NOTE:  The  position  of  zero  moment  is  a  function  of  geometry  and  the 
degree  of  fixity  at  the  ends  and  independent  of  the  type  of  loads  or 
the  magnitude  of  loads . 


LOADED  SPAN  E 
e 


Apply  three  moment  equation  to  span  D  and  E  (D  is  unloaded,  E  is 
loaded) . 


Mcd^d  +  2Mde    <*d  +  £e>    +  Mef4   +6~-4^° 

^e  \ 

Where  Re  =  Moment  of  field  load  moment  area  about  right 
support  of  span  E 
LQ  =  Moment  of  field  load  moment  area  about  left 


support  of  span  E 
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Eliminate  M  ^  and  Mp  by  using: 
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Then  we  obtain 
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from  general   equation 
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Then  above  equation  becomes: 
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Finally  we  get 
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3.   Stiffness  factor 

Stiffness  factor  is  defined  as  the  reciprocal  of  the  angle  of 
rotation  caused  by  a  unit  moment. 

From  the  elementary  mechanics  of  materials  we  know  the 
following  relation  between  angle  of  rotation  and  moments 
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By  these  relations  we  can  easily  obtain 

£  I        / 
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^CD   3EIC   rVc  6EIC 

If  we  use  the  notation  t        to  denote  the  stiffness  of  the  left 
hand  side  member,  and  substitute  ^c  by  (  ^c -  c  )  we  will  get 
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In  general,  the  stiffness  of  the  left  hand  member  at  the  right 
hand  support    £n 
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Similarly  we   can   express   in  terms   of  right  hand  member 
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In  general 


A/h  = 
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If  we  use  the  notation  Z£     including  left  and  right  hand 
members'  stiffness,  then  we  have 


A/n  = 
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Similarly  we  can  apply  the  above  formula  from  right  side  to 
left  side.   The  results  are: 
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